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Abstract 

A master equation for light generated by atoms, which 
states are prepared by a pumping mechanism that pro- 
duces atomic correlations, is derived in the Fokker - 
Planck approximation. It has been found that two- 
particle correlations only play the role under this ap- 
proach. Then the equation is applied for describing 
micromaser operations and light noise is discussed. We 
consider the correlated atomic states prepared by tele- 
cloning protocols. 

1 Introduction 

When atoms emit light statistics of radiation depends 
on the atomic noises, in particularly, a reason for which 
is correlation of atoms. To generate light atomic en- 
semble can be prepared by the various pumping mecha- 
nisms. As well known, in the case of a regular pumping 
atoms may generate light, which statistics of photons is 
sub Poisson as for the laser , and the micromaser 
operations. Thus in the case of micromaser gener- 
ation sub Poisson statistic of light can be achieved if 
atoms, are injected into cavity in such a way that only 
N atoms, where iV is a small and not random value, 
interact with a cavity mode and contributions of atoms 
are assumed to be independent. Recently it was found 
that the quasi-stationary micromaser operation by N 
independent atoms has nontrivial dynamics Q 

The main aim of this paper is to obtain master equa- 
tion for radiation, when one cavity mode interacts with 
atoms which state is prepared by a pumping that cre- 
ates correlations between atoms. Next, as example, we 
consider how statistics of generated light depends from 
atomic correlations. For this we derive a field master 
equation for the Glauber-Sudarshan quasiprobability 
in the Fokker - Planck approximation to describe in- 
teraction of one cavity mode with N two-level corre- 
lated atoms. It is not the case for using the classical 
Scully-Lamb approach [5| (see also its application in 
the micromaser theory \t\) because it starts from the 



single atom task. 

States of correlated atoms can be prepared by dif- 
ferent ways. In particularly, any pumping mechanism 
can be considered as a cloning process, which output 
is a set of N copies of an initial atomic state having an 
inversion of population of working levels. In quantum 
information processing it is already well known the pro- 
tocols of obtaining N copies of quantum states from the 
M originals. It is the transformations M N, where 
M < N. First of all this protocol can be achieved 
by a set of unitary transformations, known as Univer- 
sal Quantum Cloning Machines (UQCM) §, g, |]. 
Next, combining UQCM and teleportation, it has been 
proposed telecloning that allows a producing of 
spatially separated copies from originals by the tele- 
portation protocol. Let AI — 1, then for a pure state 
of a qubit with the help of UQCM or telecloning one 
can obtain the N - particle entangled state of the GHZ 
(Greenberger - Home - Zrilinger) class 
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where \a\^ + |/3p = 1, is the tensor product 

|6) (X) . . . (g) \b), b = 0, 1. In the same time copies of 
some quantum states can be created by teleportation 
without UQCM, when a multipaticle quantum chan- 
nel is permitted By the multiparticle channel, 
also named the multiuser channel, an initial state is 
distributed with N parties, each of which has an op- 
timal copy of the original in his hands. It is kind of 
telecloning. If the original is a mixed state, such type 
telecloning allows preparing the N - particle mixture 

Ao|0>(0|+Ai|l)(l|->Ao(|0)(0|)«^ + Ai(|l)(l|)®^ (2) 

where Aq + Ai = 1. Note that in contrast to ([l]), the 
obtained A^-particlc state (2) is separable. 

According to the no-cloning theorem all copies 
are not independent, except the case a = or Aq = 0, 
when the state of N independent particles of the form 
is obtained. Let |0), and |1) are the lower and 
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the upper atomic level, then by telecloning one obtains 
the correlated states of N atoms. The collection of 
these atoms can generate light, if they have inversion 
of the population of the working levels, when > 
|ap, Ai > Aq. Both states, given by (||) and (^, are 
quite different, but, with respect to any set of the K < 
N particles their features are similar, if iV > 2. The 
ii'-particle density matrix has the form /(12 . . . K) = 
r(|0)(0|)®^ + (l-r)(|l)(l|)®^, where r = Ao,|a|2. It 
means, that any atom has no coherence on the working 
transition, its state is mixed, and atomic correlations 
are purely classical. Next we also consider the state 

\Z) = a\bb) + (3\^+) (3) 

where \a\^ + \p\^ = I, -^1+ = (|01) + |10))/V2 is one 
of the Bell states, b = 0,1. This is entangled state, 
for which the one-particle density matrix reads /(I) = 
\a\^\b){b\ + \P\^-l/2+{aP*\b){l\Sbo + \b){0\Sbi+h.c.)/V2. 
Indeed, there is the atomic coherence, or polarization, 
that results in some new features as, for example, the 
laser and micromaser operation without inversion of 
population on the working levels fl^ . 

What is a difference between the presented corre- 
lated states, prepared by telecloning, from the states, 
produced by any "usual" pump? As for the simplest 
model of a pump, one can consider a collection of N in- 
dependent atoms, say in the state /(I . . . N) = /(l)®^, 
/(I) = Ao|0) (0 -I- Ai|l)(l|. Let atoms are illuminated 
by a strong coherent field. This interaction can be 
described by a unitary transformation, that results in 
the inversion of population of the working level. But 
one finds, that any correlation between atoms cannot 
be produced. Then atoms become active, but inde- 
pendent, its entropy has the form S = NHi, where 
Hi = — Aq log Aq — Ai log Ai. The entropy of the pure 
states (P and (^) is equal to zero, the entropy of the 
mixed state (H) is 5* = Hi. Because under an unitary 
transformation entropy does not change, the correlated 
states, we consider, can't be prepared by a simple way, 
when independent atoms interact with a classical field. 

To derive the master equation, we use formalism pre- 
sented in jl^ for considering the interaction of N corre- 
lated atoms with a quantum electro-magnetic field. In 
this approach we use the standard procedure of illumi- 
nating of atomic variables assumed to be fast. It results 
in the Fokker - Planck equation (FPE) as example for 
the Glauber - Sudarshan quasiprobability widely used 
for the quantum optics problems. In spite of the states 
(|l|), (H) and (^, we discuss here, describe entangled 



and correlated atoms, in fact it is not a real interac- 
tion between atoms and they are treated as a non in- 
teracting set of particles. So that the correlation gives 
no additional difficulties in deriving the FPE. One of 
the main result followed from the used approximation 
is that all diffusion coefficients of the obtained FPE 
involve the two-atom correlations only. Then for our 
purpose it needs consider not less then 2 atoms, in- 
teracted with the cavity modes. In the same time we 
assume the number of atoms is limited, to consider the 
scheme close to micromaser. It is in accordance with 
the standard approximation of small fluctuation of field 
that we also use. The reason is that, the diffusion coef- 
ficients of the FPE have the contribution proportional 
to N{N — 1), due from the initial multiatom problem, 
we consider here. Under the small fluctuation approxi- 
mation the diffusion coefficients must be small, that is 
the restriction on the number of atoms. 

Using the found FPE, for the case of initial state of 
atoms (3) we calculate the Mandel parameter of the 
light and its spectrum of noise for the micromaser op- 
eration. We find here that atomic correlations lead to 
increasing of the the light noise, but in particular case 
the noise can be reduced up to zero, if all atoms are pre- 
pared in the upper states Note that the perfect 
reduction of the micromaser noise can be achieved also 
by the regular pumping as it has been demonstrated in 
ref . 1^ , when the pump results in the Fock state of the 
light. However this generated state has no well defined 
phase because of its diffusion. In contrast, if atoms are 
prepared in the entangled state, given by (||), we found 
a regime in which the light is generated without inver- 
sion of population of atomic levels and in particularly, 
the generated mode has steady-state phase and noise 
suppressed up to zero. 



2 Initial equations 

Let consider a collection of N similar two level atoms 
interacting with one cavity mode having a frequency 
equaled to an atomic transition frequency. This system 
is described by the density matrix F that obeys the 
following equation 

(4) 

V = -ihg{Sioa - ^oiof) 
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where g ~ dfi~'^ {huj / {2L'^tQ) is the coupHng con- 
stant, d is a dipole moment of a working transition, 
a\ a are photon operators, [a;a^] — 1; atomic opera- 
toi S Sxy are defined as 



N 



^xy — ^ ] Sxy{k) 



(5) 



k=l 



Sxy{k) = \x)k{v\ 

x,y = 0, f 

where |0)fc is a lower and |f)fc is a upper state of an 
atom k. The relaxation processes can be taken into 
account by adding appropriate terms in (^). 

To derive a master equation for field based on the 
formalism presented in |15| , let introduce the represen- 
tation of the density matrix F over the coherent states 
of field 



(6) 



F= d^a $(a)|a)(a| 



where a\a) = a\a). The matrix <I> is the operator with 
respect to the atomic variables but a c-function with 
respect to field. When $ is averaged over the atomic 
variables, one finds the well-known Glauber-Sudarshan 
quasiprobability P(a) = SpA^- An equation for $ can 
be obtained from (Q) taking into account the correspon- 
dence ^ [i?o, ^] + 9q(D$) where field operators 
a and are presented by complex numbers a, a* as- 
sociated with the field amplitudes, and by the deriva- 
tives over these amplitudes as den da* ■ The Hamilto- 
nian t?o represented the interaction of the atoms with 
a field of the given complex amplitude a has the form 
i^o = ^('S'oiQ^* ~ c.c), the term involved the derivatives 
reads da{D^) = -gd(Sm^)/da + h.c. 
Introduce the representation for $ 



$ = / + n 



(7) 



where / is the atomic density matrix for which SpAf — 
1, where subscript A denotes spur over the atoms, 11 is 
a correlation matrix, SpA^ = 0. The field quasiproba- 
bility P = SpA^ obeys the following equation 

jrP - -gjT K^oi)^ + SpAiSoiIl)] + h.c. (8) 

where the angle brackets here and later denote the av- 
erage of the atomic operators over the density matrix 
/ as (Sxy) = Sp{Sxyf) x,y = 0, 1. It is worth paying 
attention that for / we get here the closed problem in 



which the interaction of atoms with field is described 
by the Hamiltonian t?Oj where the field is given by its 
complex amplitude only. It looks as the problem of 
evolution of the atomic ensemble in a classical field 



d_ 
di 



(9) 



Using (Q) and taking into account and (||) , one finds 
the precise equation for the correlation matrix 11. How- 
ever we wish the FPE for the field quasiprobability P. 
Therefore it needs a solution for 11 in a perturbation 
theory over derivations to a, a* up to the first deriva- 
tives only, because it gives a term to the equation for 
P (||) resulting we will have derivates of an order not 
higher then two. It results in the following equation 
for n 



^{Soi-{Soi))f + h.c 



P (10) 



To derive a closed equation equation for P usually it 
is assumed that an atomic system evolves more faster 
then field and introduced a rough time scale. Let T 
is a typical time of changing of the field and consider 
the rough derivative over T as (P(T + t)- P{t))/T = 
dP/dt. Then integrating equation (H) for P over time 
from t till T + t and taking into account the fact, that 
the field is slowly changed in the scale time of T, we 
find the desired equation for P 
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dt'^{{Soi{t'))P{t) + SpAiSoMt'j) 



da 



+h.c. (11) 



This is the master equation for field. Its coefRcients 
can be found from the the problem for / and 11 where 
the field is presented by its quasiprobability P and as- 
sumed to be constant in time. Also it needs to denote 
the initial matrices /(O) and 11(0) are taken at an ini- 
tial moment of time, say t = and it is used as a way 
to introduce the pumping mechanism. As for exam- 
ple, consider atoms placed into cavity, which Q-quality 
is modulated periodically, and atoms interact with the 
high Q cavity mode during time T. Then we can de- 
scribe a pumping as initial state of atoms, prepared at 
time t as f(t), and for an appropriate choosing of the 
time intervals we set Il{t) = 0. Also our scheme can be 
considered as micromaser, when a beam of atoms, that 
states are given by pumping, is injected into cavity, in 
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which N atoms interacts with one mode during time T 
and after leave the cavity. As result, the master equa- 
tion describes the field behavior, arisen from summing 
over the all changing of the field in time T. Note, all 
contributions are independent. 



3 The atomic averages 



According to evolution of the electromagnetic 

field depends on the atomic observables of two types. 
They are based on the atomic matrix / and the cor- 
relation matrix H, like {Sp{t)) — SpAiSpfit)) and 
SpA{Spn{t)), where p = 0, 1, 2, 3 or p = 00, 01, 10, 11, 
if the binary notation of p is used. Because atoms 
are assumed to be identical, one finds, that {Sp(t)) = 
N{sp{t;l)), and the one-particle operator /(t; 1) = 
Sp2...N f{t;l ■ ■ ■ N) is needed. The matrix / can be 
obtained by solving that looks as a problem of 
atoms in the field of the given complex amplitude a, 
in other hand it can be treated as a task of interaction 
of atoms with a "classical" field. The matrix /, as well 
as n, cannot be factorized, because of the pumping, 
we consider, prepares the correlated atoms, in contrast 
the case of independent particles, discussed in Ref [p^ , 
when both / and 11 are product of the one-particle op- 
erators. It is important, that there is no any interaction 
between atoms, in spite of the atomic correlations are 
created by the pumping. For the problem ^ it results 
in operator evolution has the factorized form and one 
finds 

/(t; l...N) = U'^^f{t = 0; 1 . . . 7V)(C/t)®^ (12) 
U = fi + i^sio(l) - iy*soi{l) 

where /i = cos{g\a\t), v — —{al\a\) sin(|a|i). 

There is another situation for calculating the observ- 
ables associated with the correlation matrix IT. In- 
tegrating (|l^), one finds that all values of the type 
SpA(Sp^{ty) depend from the atomic variances involv- 
ing pair of the atomic operators 



SpA{Sp\i{t)) = 



-.9 / dt' 

'0 



d 



da 



Dpoi{t,t') + —Diop (t',t) 



Pit) 



where the Heisenberg operators are introduced as 
Sp{t) = {W)'»^Sp C/®^ and averaging is taken over 
the initial state of atoms f{t = 0; 1 . . . N). Note, that 
in ( p^ we suggest, that the atomic operators commute 
with all derivatives over a. That is a typical approxi- 
mation for the FPE ||l^ . The variances, given by (|l^) , 
consist of two parts, where first is the one-particle con- 
tribution proportional to N, and second is the two- 
particle contribution proportional to N{N — 1) 



Dpg{t,t')^N[{sp{t;l)sg{t';l)) 
-{sp{t;l)){s,it';l))] 
+N{N~l)[{sp{t;l)sq{t';2)) 
-{sp{t;l)){s,{t';l))] 



(15) 

(16) 
(17) 



Note, in the classical Lamb-Scully approach one term 
from ( p^ N{sp{t;l)sq{t';l)) takes place. For calcu- 
lating the atomic variances, the Heisenberg picture is 
suitable, it results in evolution of the one-atom opera- 
tor in the form 



Sp(t;l) = ^i?p,(t)s,(l) 



(18) 



where g — 0, 1, 2, 3, and the unitary matrix Rpq reads 



—fj.v 
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y \v\ /ij/ ^v* 

Then all variances take the form 



(19) 



Dpq{t,t') =Y,Rpp{t)RqQ{t')DpQ 

PQ 



(20) 



where p,q,P,Q = 0,1,2,3 and Dpq = Dpq{0,0) is 
an initial correlation function, given by ([T^), in which 
t — t' — 0. It can be seen, that the two-particle density 
matrix f{t = 0; 12) is needed to find these variances. It 
means, that the two-atomic correlation only is required 
to specify the field behavior under the Fokker - Planck 
approximation. 



(13) 



where variances Dpq(t, <'), p,q = 0, 1, 2, 3 have the form 
Dpq{t,t') = {Sp{t)Sqit')) - {Sp{t)){Sqit')) (14) 



4 The Fokker-Planck equation 
and noise 

After introducing the polar coordinates / = jap, = 
arg a the master equation for field (O) takes the FPE 
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form 



dt 



d 1 
^ 9^^" +2 



U.V — I,ip 



2 dudv 



Quv 



p 



(21) 



All coefficients at the first derivatives are denoted 
by the one-particle atomic polarization {sQi(t; 1)) 



Ai = -N- 



T 



dtexp{-iip){soi{t;l)) + h.c. (22) 



2\/7t JO 



i exp(— i(/3)(soi(t; 1)) + h.c. 



At this step the field relaxation has to be taken into ac- 
count, because the light must leave the cavity. For the 
model of the Lindblad type, it results in Aj Aj + CI, 
where C <^ 1/T is a rate at which photons come out of 
the cavity. C depends on the output mirror transmit- 
ting. The diffusion coefficients are based on the atomic 
variances and involve the two-particle correlations 



Quv — 
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•t 



dt / dt'q.,,{t,t') (23) 

JO 



where 9uv = Ovu, Oi^ = 1, 9^^ = 1/(4-^). = I> 

qiip = q-fi = -i exp{-2if)Dii{t, t') -f c.c. (24) 
q^v = D2i{t, t') ~ exp(~2iLp)Dii{t, t') -f c.c. 
qii = D2i{t, t') + exp(-2iv?)Dii(t, i") -f c.c. 

To describe nose of the light we use the Mandel 
parameter ^, that indicates how statistics of photons 
deviates from the Poisson statistics: (n^) — (n)^ = 
(n)(l -|- ^), where n = a^a is the photon number op- 
erator. In the presented representation ^ = (e^)/(/), 
where e = / — (/) is a fluctuation, that is the differ- 
ence between "intensity" of the light / and its aver- 
age value (J) = J d'^P{a)\a\^ . Assuming, that fluc- 
tuation of intensity is small and does not depend on 
the phase fluctuation, then it follows from (|l]) imme- 
diately, that the Mandel parameter is proportional to 
the diffusion coefficient Qn as ^ = Qii{{I)T)^^ , where 
T = C + [dAi / dl)i^(^i) is the decay rate of the in- 
tensity fluctuations. The considered approximation of 
small fluctuation is the standard one for analysis the 
FPE and it needs ^ < (/> • 



5 Noise of micromaser operation 

To apply the obtained FPE consider the micromaser 
operation, when the atomic relaxation can be ne- 



glected. It is true, if the decay rates of atomic lev- 
els 7 such that 1/T » 7 C. Assume fluctuations 
of intensity are small in a sense, that the deviation 
e = /—(/) <C (/) slightly change near a steady state 
solution (/). 

First consider the pump, produced classically corre- 
lated atoms in the state, given by (|l|), if iV > 2 and 
(1). Then from (|2^) it follows, that (/) is the steady 
state of the semiclassical micromaser equation of the 
form 



|(/>^-C(/> + ^(Ai-Ao)sin2i? 



(25) 



where B = g^{I)T. For steady state it needs an 
inversion of population Ai > Aq and the condition 
(1 — S/tanS) > due from the requirement of the 
stable equilibrium. 

Then using ( |23| ) under small fluctuations of intensity, 
one finds the Mandel parameter 



T 



Q 



II 



N (2Ai - l)sin^S(l - B/l&nB) 



where the diffusion coefficient Qn has the form 



Q 



II 



(26) 



(27) 



N 1 

— {--sin^B-l- (2Ai - l)BsinScosB 



(1 - Ai)sin^B 



+2Ai(l-Ai) sink's 
-2(iV- l)Ai(l - Ai)sin'^B} 



In ( p7| ) the atomic correlations of both state (|l|) and 
(H) presented by the same term N2{N — l)Ai(l — 
Ai)sin'^fe/r are non negative. It means increasing of 
the light noise due from the initial atomic correlations. 
This term is equal to zero, if = 1, or, if the pump 
creates independent atoms. Indeed, the phase behav- 
ior is insensitive to these atomic correlations. In ac- 
cordance with the small fluctuation approximation it 
needs Qjj ^ r(/)^, that is the reason of restriction 
the number of atom N. 

Consider the noise of light, denoted as the spectrum 
of photocurrent (w) . It can be measured by a simple 
heterodyne scheme included a detector and spectrum 
analyzer. At frequencies w « the spectrum of the 
light noise takes the form 



,■(2) 



(0) = 1 + 2e 



1 - B/tanB 



(28) 
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where unit indicates the shot noise level or the standard 
quantum limit, which is noise of the coherent field. The 
main features of the light statistics can be found from 
equation (p6[),(|2^). If all atoms are prepared indepen- 
dently in the state then ^ > -1. The negative 
Mandel parameter means, that in the cavity the pho- 
ton statistics is sub Poisson, and the noise of light can 
be suppressed below the standard quantum limit. It 
^ = — 1, the photon number variance is equal to zero, 
and one finds the Fock state of the light into cavity \n), 
where n = (/). However the noise reduction is small. 
Noise can be suppressed up to zero, when — —1/2. 
These results are in agreement with results for a reg- 
ular pumping, proposed in Ref. |^. For the case of 
correlated atoms, we discuss here, the Mandel param- 
eter becomes positive and the light noise increases. 

In contrast, when a pump prepares atoms in 
the state given by (||), the atoms have the ini- 
tial coherence (soi(O)) = {aP*6bi + a*/3Sbo)/V2 = 
|a/3| exp[?(po], where ipo = (— l)^^"**^ arg(a/3*) and in- 
version (sii(O)) - (soo(O)) i-iy-''\a\^. Here b = 
0, 1, that indicates two version of the state (|). If 6 = 
the inversion of population is negative, in contrast the 
case when 6=1. Then together with the small fiuctu- 
ation of the intensity, assume the phase fiuctuation is 
small in a sense, that n = tp — ipo <^ ipQ. Then semi- 
classical equations for the phase ipo and the intensity 
(/) take the form 



d_ 

dt 



ipo = -V2-^\a(3\ sin('0o - (Po) 



|(/) = -C(/) + |;(-l)i-Vl'sin^i? 
2 

+ — cos(t/'o — ^o) sin 5 cos B 



(29) 
(30) 



In the steady state V'o = fo ^-nd this solution is sta- 
ble. The initial atomic coherence is presented by the 
last term in the right hand side of (|3^). It result in 
the steady state solution for intensity, even the initial 
inversion of population, produced by the pump, is neg- 
ative, when 6 = 0. 

Noise of the light or spectrum of the photocurrent 
at the frequencies near zero can be written as 



z2(0) = l + 2e^ 



(31) 



where 



Y^C(1~ s ('^('t9B + V2\aP\{ctg^B-l)) 
{A + V2\a(3\ctgB) 



where A = {-l)^^-'>h\a\'^ , ^ = Q///((/)r). The ex- 
pression for the diffusion coefficients reads, if 6 = 1 

Qji = ^[+^A{sm'^ B + 2BsmBcosB) (32) 



if 6 = 
Qii 

and 



\/2|a/3| cos 2B{B - sin B cos B) - Q] 



= ^[^Aain^ Bisin^ B + 2)-2BcosB 

+ V2\a(3\cos2B{B + smBcosB) 

- 2 sin B cos B - Q] (33) 



2 sin^ B cos^ B (34) 



Q = (2(H-2|a|' 

+ V2\a/3\AsmBcosB{l - cos2B) 

- -A' sin^ B 
2 

There are two types of steady states of the micromaser 
operation, that can be pointed. For the first of them 
6 = 0, it means, that atoms are prepared initially with- 
out inversion because (sii(O)) — (soo(O)) = — If 
\a/3\ « 0.49, then one finds light, which phase is well de- 
fined and statistics of photons is sub Poisson. Indeed, 
the light noise can be suppressed up to i^{0) w 0.27, 
that is limit here. The second type of operation is ob- 
tained, when 6=1. In this case the pump creates 
inversion of population of the atomic levels and co- 
herence. Atoms generate light, which level of noise is 
various. Particularly, one finds the perfect noise re- 
duction, when i^{0) ~ 0, if jap '3> |/3p. Physically, the 
last feature is clear, because the initial state of atoms 
is close to a state, of two independent atoms in the 
upper level. Indeed, in the micromaser operation the 
ideal sub Poisson light can be achieved using a regu- 
lar pumping of atoms However note that to do it, 
atoms have to be pumped by the light, which statistics 
of photons is already sub Poisson ||l^ . 

This work was supported in part by the Delzell Foun- 
dation Inc. and RBRF grant committee. 
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